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I We systematically explore the parameter space of the state-of-the-art brane- 

I antibrane inflation model (Baumann et al, arXiv:0706.0360, arXiv:0705.3837) which 

, is one of the most rigorously derived from string theory, applying the cosmic back- 

. ground explorer normalization and constraint on the spectral index. We improve 

on previous treatments of uplifting by antibranes and show that the contributions 
CN ■ from noninflationary throats play an important role in achieving a flat inflationary 

potential. To quantify the degree of fine-tuning needed by the model, we define an 
effective volume in the part of parameter space which is consistent with experimental 
^, constraints, and using Monte Carlo methods to search for a set of optimal parame- 

^ . ters, we show that the degree of fine-tuning is alleviated by 8 orders of magnitude 

relative to a fiducial point which has previously been considered. In fact, close to the 
CN ■ optimal parameter values, fine-tuning is no longer needed for any of the parameters, 

^r^ I We show that in this natural region of the parameter space, larger values of Ug close 

I to 0.99 (still within 2a of the WMAP5 central value) are favored, giving a new aspect 

of testability to the model. 

o ■ 

O '. I- INTRODUCTION 
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Brane-antibrane inflation is one of the most distinctive and highly developed applications 



^ [ of string theory to cosmology. Although the basic idea seems simple — getting inflation from 
the potential between a brane and antibrane which are separated along an extra dimension — 
making it work has proved to be quite challenging. The earliest versions were incomplete due 
to the lack of an explicit mechanism for stabilizing the moduli, notably those associated with 
ihe size and shape of the extra dimensions. Major progress in this respect was made in Ref. 
l|, where warped compactification of type IIB string theory using the Klebanov-Strassler 
throat [2I was exploited. 

In Ref. [l| , it was shown that the potential for a brane falling down the throat suffers from 
the rj problem: the curvature of the inflaton potential is generically of order H^, making it 
too steep for inflation. To counteract this, the generically large inflaton mass must be tuned 
to a small value using a canceling contribution from the F-term potential, coming from the 
superpotential 

W = Wo + Ao(0)e-'^^, (1) 
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where (f) denotes the brane position (inflaton) and T is the complex Kahler modulus, whose 
real part a determines the overall volume of the extra dimensions. The second term arises 
from nonperturbative physics like gaugino condensation on a D7-brane, which is one of the 
essential ingredients of the compactification. The existence of the F-term potential was 
already known to be necessary for stabilizing T [s^; however the dependence of Aq on was 
not known by the authors of Ref. they merely parametrized this dependence. 

The actual dependence Ai((/>) however is calculable within string theory, and this missing 
step was carried out in Ref. ^. It was subsequently shown ^ that the explicit form of Ao(0) 
was apparently not amenable to achieving the desired cancellation to get a small inflaton 
mass. However this conclusion depended on exactly how the D7-brane was embedded in the 
throat. A different embedding was considered in Refs. gI, Q, S] which yielded the desired 
form of the F-term potential. With sufficient tuning of parameters, it was possible to achieve 
inflation. Because of this need for substantial tuning, the scenario was dubbed "delicate" in 
Refs. 0,0]. The difficulties were further elaborated in Ref. which surprisingly concluded 
that it was difficult to satisfy both the cosmic background explorer (COBE) normalization^ 
and the constraint on the spectral index. Normally one expects that there is enough freedom 
to adjust the overall scale of the potential to achieve the normalization of the amplitude, 
and indeed, we will show that satisfying the normalization presents no extraordinary new 
difficulty.^ 

On the other hand, we do point out an additional challenge which was not consistently 
treated in previous work. Namely, it is necessary for the cosmological constant to vanish at 
the e„d of inflation. In Ref. i, sufiieient flatness of the potential was achieved by tuning the 
parameter s, which is the ratio of the uplifting energy to the absolute value of the negative 
energy of the anti-de Sitter (AdS) minimum in the absence of uplifting. In fact, this ratio is 
fixed (to approximately unity) by the necessity to uplift to a Minkowski minimum, thereby 
removing it as a tunable parameter. We take advantage of the fact that its tunability can 
be restored by considering the contributions from other throats to the uplifting. 

In Refs. [tI, [Sj, it is claimed that the brane-antibrane inflationary model must be rather 
finely tuned. However, there is little quantitative basis for this statement. Part of our 
objective is to invent a quantitative measure of the degree of tuning. We propose a measure 
of tuning via the relative volume of parameter space around a given point, which is consistent 
with inflation. This statistic allows us to say to what fractional precision the parameters 
have to be adjusted in order to get inflation consistent with observations; e.g., parameter x 
must be tuned to one part in 100, parameter y to one part in 1000, etc. Using this measure, 
we show that indeed, the model is quite flne-tuned at the parameter values which have 



^ The historical term "COBE normaUzation" refers to the value of the primordial power spectrum at 
k 0.002 Mpc^^ of course, we use the latest Wilkinson Microwave Anisotropy Probe (WMAP) value. 

^ We will see in Section IV Bl that the normalization and the spectral index are correlated; hence further 
fine-tuning is needed to construct a model satisfying the current observational constraints. 
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been previously suggested. However, we will show that there exists a more favorable set of 
parameters around which the tuning problem is much less severe, by adapting Monte Carlo 
methods which have been widely used in other cosmological applications for this purpose. We 
also consider the degree of fine-tuning using an alternative measure — the absolute volume of 
the allowed parameter space; we compare the two measures and show that the improvements 
in fine-tuning are of the same order. 

A further delicate aspect of the model is the problem of overshooting the flat part of the 
potential by starting too high 1^. With the wrong initial conditions, the infiaton gathers 
too much speed to roll slowly in the fiat region of the potential. We investigate the scope of 
initial conditions which is compatible with inflation and show that this aspect of the tuning 
problem is also improved somewhat in the most favorable parameter range. 

A final technical issue concerns the extent to which the model can be treated in terms of 
a single field, due to the fact that the trajectory curves in the space of 0, the brane modulus, 
and T, the Kahler modulus. A seemingly reasonable treatment, involving both fields plus an 
analytic approximation for the dependence of T on 0, can result in an incorrect prediction 
for the spectral index. We show that the single-field approximation can be a good one, 
provided one is sufficiently careful. 



II. THE POTENTIAL OF L>-BRANE INFLATION 

The infiaton potential derived by Baumann et al. 0, S] consists of two contributions: a 
supersymmetric (SUSY) F-term potential Vp-, and a SUSY-breaking contribution due to 
anti-D3-branes, Vd, which is necessary for uplifting the minimum from a negative value (an 
AdS minimum) to zero. Vp depends on the correction to the nonperturbative superpotential, 

/lo(0) = Ao^?^/"(0)=Ao 

Aq is now taken to be a constant, n is the number of D7-branes in the stack, and cf)^ is 
proportional to the radial position where the D7-brane stack comes closest to the bottom of 
the throat. For < 0^, the D3-brane is below the D7-brane stack in the throat; this is the 
region where inflation takes place. The configuration is illustrated in Fig. [H 

The dynamical degrees of freedom are expressed in terms of the dimensionless rescalings 
of the D3-brane modulus and the Kahler modulus cr. 




oj = aa = — a. (4) 
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FIG. 1: Schematic representation of D7-brane stack and D3-brane in Klebanov-Strassler throat. 
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where c = 9 / (AnujQCp'i) , and Uq is defined to be the stable value of cu when a; = 0, i.e., 



dV 

duj 



= 0. (8) 

Thus ujq is not a free parameter. (See Appendix lA II for details.) The functions U and D 
come from the Kahler potential and the brane-antibrane Coulombic interaction, respectively, 
and are given by 



U (x, uj) 
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(9) 
(10) 
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The uplifting term D(x) is generalized relative to Ref. [8|] by including an extra contribution 
Di representing the contribution to the uplifting which remains after the brane-antibrane 
annihilation.^ The first term involving the warped antibrane tension Dq in the infiationary 



^ We set the reduced Planck mass to Mpi = 1. 
The revised versions of Ref. [8| have -Dothcr in Eq. (2.20) which plays the role of Di. However, the authors 
still vary the total uplifting to tune their potential, rather than the ratio of the upliftings. Comparing 
Eq. (3.19) in Ref. ^] to Eq. p^. you will see that we allow the ratio of the upliftings to change. 
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throat has been resummed rather than Taylor expanded in the Coulomb interaction. The 
difference between these two forms is negligible in the slow-roll region of the potential, but 
the resummed form has good behavior as x — 0, unlike the Taylor-expanded version, and 
this is convenient for numerical evolution going to the end of inflation. (In the region where 
the potential is fiat, the Coulomb interaction typically plays an unimportant role, and it is 
sufficient to approximate D ~ Dq + Di.) The resummed form resolves an inconsistency in 
the treatment of Ref. [sl, which relies upon the value of the potential at a; = 0; however this 
only exists in the expanded form if one ignores the Coulombic part altogether. 

Since the x — >■ limit of the uplifting term is not treated consistently in most of the 
existing literature, it is worthwhile to give some further explanation.^ In a more realistic 
description, we would need to include the tachyon field of the brane-antibrane system, whose 
mass squared becomes negative when the separation between the two falls below some critical 
value of order of the string length scale. Instead of x ^ 0, the evolution would continue in the 
tachyonic direction. For our purposes, these details are not important because inflation has 
already ended by this time. What is important however is to know how much the potential 
decreases between inflation and the minimum of the potential, and the latter must be at 

= so that there is no residual cosmological constant. Our resummed expression, which 
goes to zero as x ^ 0, correctly models the fact that the tension of the brane-antibrane 
system is precisely the amount by which V changes between inflation and the minimum 



12l |. Although Ref. 12] worked in an unwarped background, it is obvious that even in 



the warped case, the uplifting provided by an antibrane must precisely vanish once it is 
annihilated by a corresponding brane. This argument shows that our new parameter Di 
represents the residual tension of the branes left over after the annihilation, assuming they 
are in different throats from that of the inflationary brane. If they were in the same throat, 
Di would have to be an integer multiple of Dq, counting the number of antibranes in the 
stack which remains after annihilation of the mobile brane; this would reduce our ability to 
tune parameters. [Also Cd in Eq. ffTOl) would have the extra factor 1 + Di/ Dq due to the 
Coulombic attraction between the mobile brane and the antibranes in the stack.] On the 
other hand, antibranes in separate throats have a tunable tension, via the warp factors of 
the throats, since it is the warped tension which appears in the potential.^ 

In addition to the D3-brane's radial position in the throat, there are five angular direc- 
tions, as well as the imaginary (axionic) component of T. Some of these have large masses 
and have been set to the values which minimize their potential. Others have a nearly flat 
potential due to approximate isometrics of the throat geometry, but as usual, such compact 
directions cannot give rise to any significant amount of inflation because their motion is 
quickly Hubble damped. (For interesting possible effects of these fields on the generation of 



^ See the appendix of Ref. ^] for a derivation. 

^ Of course the warp factors themselves are determined by exponentiated ratios of integers; we assume these 
can be approximated by continuously varying warp factors. 
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density perturbations at the end of inflation, see Ref. ISj.) On the other hand, the Kahler 
modulus uj can undergo significant evolution, so it is important to keep both the uj and x 
fields at the outset, even though only one linear combination is light at any point along the 
infiationary trajectory. 

We see that the potential depends on six parameters, Aq, 0^, n, Wq, Dq, and Di. To 
make closer contact with Refs. p, 0], we will use a different parametrization in place of Wq 
and Di. The stable value of the Kahler modulus before uplifting, up, is defined through 



which gives 
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Therefore Wq can be traded for up- In place of the parameter Di, we define the ratio of Vd 
and Vp- 
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(13) 
(14) 

(15) 
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where Dqi is the ratio of two tensions: 
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As a result, the potential is determined by the following six independent parameters:' 

s,Ao,Doi,(f)^,ujF,n. (18) 

The dependent parameters a, c, and Cd have been defined above. 

These parameters of the potential can be expressed in terms of more fundamental string 
theoretic quantities which imply constraints on their allowed values. In particular, 0^ 



^ In Section [V CI we will show that s is a function oi ujp when the requirement of uplifting the potential is 
imposed. 
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and uip can be expressed through stringy parameters: 

2 



(19) 
(20) 



where > 1 is the 5- form flux (equal to the product of the 3-form fluxes), > 1 is 
the ratio of the bulk part of the 6D volume to the part due to the throat, i?4 > 1 is the 
corresponding ratio for the 4-cycle volume on which the D7-brane is wrapped, and > 1 is 
the ratio of the radial coordinate ruv at the top of the throat to the value corresponding 
to 0^ in Fig. [H (See Appendix IA2I for the complete list of theoretical constraints on the 
microscopic parameters.) We will take these constraints into account when we search the 
parameter space. 

In addition to the above, there is another important constraint which was not discussed 
by previous papers. The parameter Aq is related to the scale A of gaugino condensation by 
Aq = A^ (see, for example, Ref. ^). This should certainly be less than the Planck scale, so 
one should at least demand that Aq < 1, and more realistically it should be even smaller. 
Below we will find a preferred value of Aq ^ 0.01. 

The shape of the potential is shown in Fig. [21 for the optimal parameter set ( l60l) which 
will be discussed later. Away from the inflection point near x = 0.03, it becomes steep, and 
starting from very large values of x can lead to an overshoot of the flat region. 
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FIG. 2: The potential for the optimal parameter set ()60p (to be discussed later), in the vicinity of 
the inflection point. The inset shows the potential over a larger range of x. 
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III. EQUATIONS OF MOTION 

The kinetic term of the D3-brane comes from the Dirac-Born-Infeld (DBI) action, of the 



form —a^T^yl — a'^X^, where a is the warp factor in the throat. For small enough field 
velocities, we can expand to leading order in to obtain a conventional kinetic term. We 
will check the validity of the approximation a~^X^ ^ 1 in our numerical analysis. In the 
framework of the low-energy supergravity effective action, the kinetic term for the D3-brane 
and the Kahler modulus comes with a nontrivial metric Gij on the space of the complex 
fields, taking the form: 

>Ckin = Grjd^^'d^^'. (21) 

This leads to the kinetic term: 

Ain = ^(3c^' + 2ujQUj(t>lx'). (22) 
We define the canonical momenta: 



dC AujQLJ(f) 



X, (23) 



duj a^U^ 
which can be solved for the field velocities. 



dx a?U'^ 
9^ 6 . 



= —^T^co- (26) 
D 

The equations of motion for the fields are 

d 

TTi + SHiii = —— (£kin -V) , i^x,uj, (27) 

d(/)i 

3H^ = Cun + V. (28) 
It is convenient to use the number of e-foldings instead of time as the independent variable: 

^ = H^. (29) 
dt dN ^ ' 

Then the equations of motion, in suitable form for numerical integration, are 
dN = -'"^ + ^i^^^^^^-a^J' 

- -37r, + - -— /:kin+ o,,o - ^ h (31) 



dN H\ aU '^'^ a?V du 

dx 
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where x and u are determined through Eqs. (125|) and (l26ll : the derivatives of the potential 
are given in Appendix IA3[ 

In the regime that we are interested in, near the tip of the throat, 3> a; and a; ~ tuo; 
then the equations of motion can be approximately written in canonical form 

^kin = + \x\ (34) 

+ 3/f0 = -V^, (35) 

X + ^Hx = -Fx, (36) 

where 




and Vx denotes dV/dX. Although we will solve the exact equations of motion in this 
paper, the canonical fields are useful when discussing the primordial power spectrum and 
the spectral index, as will be seen in the following sections. 



IV. THE PRIMORDIAL POWER SPECTRUM 



The primordial power spectrum for single-field inflation is 



Vn = (38) 



Generalizing to the noncanonical kinetic term, it can be written as [l^] 
The spectral index is 



= (39) 



-1 



dlnVn d\nH dlnC^^^A d\nH . 

"^"' + ^h^ = '+i'^iv rfi^Ji' + ^j ' 

where dH/dN and dC^i-a/dN can be calculated through the potential and its derivatives 
(see Appendix lA 41 for details). Equation (l39ll is justified when there is effectively just a 
single field contributing (single direction in field space), since it is invariant under field 
redefinitions. However, it is not valid when the entropy perturbations make a significant 
contribution to the power spectrum. 

The power spectrum and the spectral index can be expressed in terms of the slow-roll 
parameters by defining the adiabatic direction if) in field space, tangent to the infiaton 



trajectory [l5| 



-0 = 0cos6' + xsin6', (41) 
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where (j) and x assumed to be canonically normalized, and 



cos^ = (42) 



sine = I ^ (43) 

li 9 7^ 0, then the trajectory is curved and entropy perturbations can source curvature 
perturbations on large scales (A; — ^ 0). However if 9 is very small, or if the entropy mode is 
suppressed by large curvature of the potential in the direction orthogonal to ip, which is the 
case in the present model [ol, the entropy mode can be ignored and then we have the usual 
formulas in terms of if), 

Us = 1 - 6e^ + 2r]^^. (45) 

The definitions of the slow-roll parameters are given in Appendix lA 51 We will see that in 
the regime we are interested in, the slow-roll approximation is well satisfied, hence Eqs. ( 139|1 
and (jS]) agree with each other. 



V. TUNING OF PARAMETERS 



Using numerical integration and Monte Carlo techniques, we have undertaken a system- 
atic study of the inflationary dynamics over the full parameter space of the model. In the 
following sections, we will first reproduce the known result jsl, that, while keeping other 
parameters fixed, the tuning of the uplifting parameter s allows one to obtain inflation with a 
sufficient number of e- foldings, > 50. We then show that by varying the amplitude of the 
nonperturbative superpotential, Aq, one can satisfy both the COBE normalization and the 
WMAP constraint on the spectral index. (Recall that Ref. claimed that it was difficult 
to satisfy both.) We next point out that the s parameter is actually already fixed by the 
requirements of uplifting; however, there remains sufficient freedom to get a flat potential 
by varying Dqi and Up- 



A. Varying the tension to get flat potential 

References jsl, ^ showed that by varying the value of Dq (proportional to the warped 
D3-brane tension at the tip) one can obtain sufficiently many e-foldings of inflation. Here 
we further investigate the dependence of the potential on s, which is related to Di through 
Eq. ( |T^ . and fix other parameters as given in Ref. j^: 

Ao = l, n = 8, Bi = 9.15, Bq = 1.5, Q^, = 1.2, = 32. (46) 
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These values imply that 



0.2406, Up = 10.009. 



(47) 



We set our new parameter Dqi to 1 for the moment for definiteness.^ We use initial condi- 
tions: 

Xi = 0.8, Ui = uj^{xi), Xi = uji = 0, (48) 
where a;* is the instantaneous minimum which satisfies 



dV 



0. 



(49) 



The initial value Xi = 0.8 is sufficient for getting inflation while avoiding the overshoot 
problem which we will discuss in Section IVll CI Taking the initial velocities to be zero is 
justified since any nonzero values would be quickly Hubble damped. 

We now consider how the inflationary potential and the resulting solutions depend on s. 
Figure [3] shows the total number of e-foldings as a function of s. The parameter space of s 
can be divided into six regions. They are^ 
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FIG. 3: The total number of e-foldings versus the value of s. Dashed lines and numbers identify 
the regions described in the text. 



This is the default parameter set in this and the following two subsections; however, we allow to vary 
in Section IVBl and {Ao,Doi) to vary in Section IVCl 

In the flat region of the potential where inflation takes place, the Coulomb term is typically unimportant; 
therefore the correspondence between D{x) used by previous authors and us is approximately D{x) —^ 
Di{l + Dqi); hence we would take Di to be smaller by a factor of (1 + Dqi) to reproduce their results. 
Moreover, since the second term gives a vanishing contribution to uplifting at x = 0, the correspondence 
between previous authors' value of s and ours is s — > s/(l + Dqi). This explains why we need s near 0.5 
for getting a flat potential, while Refs. 0, Q had s 1 for the same parameters. 
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1. s > 3.443: dV/duj ^ 0; there is no valley- shaped potential and no stable trajectory. 

2. 0.5338 < s < 3.443: the inflaton gets trapped in a local minimum, N —* oo. 

3. 0.5186 < s < 0.5338: the potential is monotonic; for s — >• 0.5186 or 0.5338, it is 
extremely flat and gives a large number of e-foldings. 

4. 0.5062 < s < 0.5186: the inflaton is again trapped in a local minimum, so ^ oo. 

5. 0.5038 < s < 0.5062: there is a local minimum in the potential; however the inflaton 
has enough momentum to escape from it if slow roll had been attained earlier (e.g., 
starting with zero initial velocities Xj = tjj = at Xj = 0.8). 

6. < s < 0.5038: the potential is negative at < a; < 0.8, so inflation ends in a big 
crunch. 

We are interested in the cases where inflation will end (although not via a big crunch). 
Although region 5 allows for the inflaton to escape from the local minimum and end inflation. 
Fig. [3] shows that one needs to extremely fine-tune the value of s to get enough e-foldings. 
Therefore we will focus on region 3, where the potential is monotonic, and much less tuning 
is needed. It is striking that there is a plateau in this region, signifying a minimum number 
of e-foldings of about 30. It is tempting to look for other parameters such that the minimum 
number of e-foldings of the plateau would be increased to more than 50. If that were possible, 
then one would say that inflation is generic rather than fine-tuned, at least with respect to 
the parameter s. To investigate this possibility, we need to search the multidimensional 
parameter space. This will be discussed in Section |VT1 

Region 3 is also the case that Ref. focused on. The latter found that by adjusting the 
tension, one can have inflation with a correct spectral index and enough e-foldings, but they 
did not succeed in finding a model which satisfies the COBE normalization simultaneously. 
We confirm this: by restricting to the default parameters in Ref. jsl and just varying the 
value of s, one cannot achieve Ug ~ 0.96 and V-ji ~ 2.4 x 10~^ at > 50. 



B. Varying A^: COBE normalization 

Of course, by varying just one parameter, one should not expect to satisfy several exper- 
imental constraints. One usually realizes the COBE normalization by adjusting the overall 
scale of the potential. In Ref. [sl, it was assumed that the prefactor Aq played the role of the 
overall scale. However this is only true as long as the Coulombic interaction is negligible, 
because as shown in Eq. (IT^ . Aq is related to Di and hence contributes to the shape of the 
Coulomb interaction through Eq. (fT6ll .^° Therefore as long as Dqi > 0, which must be the 



Even if we Taylor expand the Coulomb term as is done in Refs. [7|, |3, l9|, the shape of the Coulomb 
interaction is still related to Aq. 
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case, Aq does not merely determine the overall scale of the potential. 

This can be seen explicitly by calculating the total number of e-foldings of inflation 
versus s for varying values of Aq. If the shape of the potential did not depend on Aq, then 
the number of e-foldings would be completely insensitive to Aq, as long as the slow- roll 
approximation is valid. Figure H] shows that in fact A^tot depends rather strongly on Aq: for 
smaller values of Aq (closer to those needed for the COBE normalization), the range of s 
values which give a monotonic potential becomes smaller. An interesting by-product is that 
the minimum number of e-foldings becomes bigger. This shows that when Aq is adjusted 
to satisfy the COBE normalization, s does not need to be fine-tuned to a special value to 
ensure that the number of e-foldings of inflation is sufficient. However, s still needs to be 
within a narrow range to avoid the problem of the potential developing a local minimum in 
which the inflaton gets stuck, with no end to inflation. For example, taking Aq = 0.005 and 
s = 0.5253, we have Vv, = 2.41 x lO'^ at = 73 (A^tot ^ 6000), where = 0.944, within 2a 
of WMAP5's result ilC|.^^ We have an existence proof that it is possible to nearly satisfy all 
the cosmological constraints in the D3-D3 inflation model; however this particular example 
appears to be flne-tuned. Below we will show that less flne-tuned examples can be found. 




IQl 1 , 1 , 1 , 1 

0.52 0.525 0.53 0.535 

s 

FIG. 4: The total number of e-foldings versus the value of s with different values of ^o- 



C. Accounting for the uplifting constraint 

Up to now we have regarded s as a free parameter, as was done in the previous literature 
0, 0]. But as we have mentioned, s should already be determined by the requirement of 
having vanishing cosmological constant at the end of inflation. In this section, we take this 



" The marginalized values (mean and 95% C.L.) from WMAP5 are Vn = (2.41 ± 0.22) x 10"^ and Us = 
0.963 ±0.028. 
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requirement into account and show that it is nevertheless still possible to tune to obtain a 
fiat potential using the extra parameter Dqi, which in fact has a similar qualitative effect to 
varying s. 

First we show how the value of s is fixed by setting 

l^(0,a;o) = 0. (50) 

Using Eq. I\A2L luq = ujo{s,ujf), we find that 

{2iUF + 3)e'^°"'^^ - 2CU0 - 5 = 0. (51) 

This shows that uq can be expressed as a function of up only, uq = u!q{ujf). From Eqs. 
and fl5T]) . we have 

(52) 
(53) 



Up \ 2uJo + 5 

Up 



1+3-^ ' 

UJF UJo 



where the latter form assumes uo, ojp ^ 1-^^ Therefore, once ujp is given, both ojq and s 
are fixed by uplifting V{^,ujq) = 0, and our use of s to flatten the potential in the previous 
sections is seen to be invalid. For example, taking the value of in Eq. (H71) . we obtain 
too = 10.100, s = 1.0087. 

To compensate for not being able to vary s, we can adjust the ratio of the tensions in 
the inflationary versus the other throats, i.e., Dqi, while keeping s fixed. We find that the 
quahtative dependence of the shape of the potential on Dqi is similar to the dependence 
on s. There are five regions of the Dqi parameter space which correspond to the first five 
enumerated for s in Section IV A[ The sixth region, where the potential became negative, 
no longer exists because we have now adjusted s to avoid this problem. Furthermore, the 
dependence of the shape of the curves A^tot(-Doi) on the parameter Aq is just like that for 
iVtot(s), as can be seen in Fig. [51 Again, as Aq decreases, the minimum number of e-foldings 
of the plateau increases, while the monotonic region shrinks and more fine-tuning is needed. 

We find that varying only (-Dqi, ^o) while keeping other parameters fixed is not sufficient 
for satisfying all the experimental constraints, but varying {Dqi,ujp) allows us to nearly do 
so. For example, taking top = 15 and Dqi = 0.1227 gives V-jz = 2.41 x 10^^ and Ug = 0.949 
at = 78, out of the total number of e-foldings A^ot — 10'^. Generically we might have 
found that by varying two parameters, the two constraints V-ji ~ 2.4 x 10~^ and ~ 0.96 
could be satisfied at the COBE scale, which we take to be 50 < A^e ^ 60. The fact that 
we cannot do so here is just due to unfortuitous values for some of the other parameters 
which are not varied. In the next section we will consider variations in the full parameter 



In most cases, tuo — ^f, so s ^ 1. 
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FIG. 5: The total number of e-foldings versus the value of Dqi with different values of Aq. 



space, but here we are primarily concerned with the degree of fine-tuning needed to get 
sufficient inflation. The present example requires fine-tuning at the level of 0.13% for Dqi. 
the potential is monotonic from Dqi = 0.122 62 to 0.122 78, so the relative fine-tuning is 
1.6 X 10~^/0.1227, approximately 1 part in 1000. Below, we will identify other regions of 
parameter space where this problem is significantly alleviated. 



VI. SOLVING THE FINE-TUNING PROBLEM 



Although the required values in the previous section may subjectively appear to be rather 
finely tuned, this is a notion which requires definition. One must distinguish fine-tuning from 
the more mundane necessity of fixing parameters to agree with experimental measurements. 
In this section we propose a measure on the volume of the experimentally allowed part of 
parameter space which will allow us to quantify the degree of fine-tuning needed in any 
localized region of the space. We calculate this statistic while doing a systematic search of 
the parameters. The results are described in the final part of this section, demonstrating 
that the fine-tuning problem is ameliorated for optimal values of the parameters. 

Our goal now is to scan the parameter space in search of regions where less tuning is 
required. To make this quantitative, we need some specific measure of the degree of tuning, 
which varies locally in the parameter space. Suppose we have identified a set of parameters 
pf''^ (where i runs over the number of parameters) which satisfy some necessary criteria for 
inflation. We can then vary each parameter to find the maximum range pf'^ — Oi < pi < 
pf"* + (7j for which these criteria are still satisfied. The width of this interval is 20", . 



In general the interval will not be symmetric; it has the form p\ < Pi < p\ <jf ■ To simplify 
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A first guess for a quantity wliicli is anticorrelated witli tlie degree of fine-tuning would 
be tlie volume in parameter space consistent with inflation, defined by the product of all 
the intervals 2(Tj. The exact A^-dimensional volume would be more complicated than our 
rectilinear approximation; we ignore this, and treat the CTj's as independent quantities, i.e., 
when determining cTj, we fix other parameters at their central values, pf^ {i ^ j). The 
volume is thus given by 

a'^n-^.- (54) 



However it would be naive to think that maximizing Eq. (IMI) corresponds to minimizing the 
tuning of parameters, because this would artificially reward parameter values that happen 
to be large in absolute terms. We are really interested in the relative variation of a given 
parameter. Therefore a better measure of naturalness is the relative volume, 

'^^ = n^ = :^- (55) 

I Pi 11 Pi 

i 

We can furthermore define a reduced relative volume, 

5'^ = (56) 

whose reciprocal quantifies the average degree of tuning per parameter, and thus makes it 
meaningful to compare searches in which different numbers of parameters are varied. 



A. Description of algorithm 

To explore the parameter space, we used the Metropolis algorithm, which looks for a 



function's minimum by the method of simulated annealing 17|. In our case, we chose the 
objective function /obj to be the negative of the relative volume, Eq. fl55l) . Because we are 
interested in the part of parameter space corresponding to Figs. [3] and HI where there is 
always a minimum number of e-foldings of inflation, we chose as our criterion for successful 
inflation that the potential should be monotonic, rather than having a local minimum where 
inflation would never end. In addition, we want to select models that satisfy the experimental 
constraints. We do this by requiring that the central values pf'^ of the regions correspond 
to experimentally allowed models, although we do not impose this additional requirement 
on the neighboring points that define the volume. Specifically, for the central point of each 
allowed volume we demand that the COBE normalization {V-ji = 2.41 x 10^^) is satisfied 



computations, we take ai to be the minimum of these two values, which underestimates the allowed 
volume. 

^'^ On the other hand, it can also be argued that the absolute volume is also a reasonable measure of the 
tuning, and we will come back to consider it in Section IVl CI there we will show that our conclusions are 
not sensitive to this choice. 
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with the number of e-foldings 50 < < 60 before the end of inflation. We take this to be a 
reasonable reflection of uncertainties in the time of horizon crossing due to variations in the 
scale of inflation and reheat temperature. We also require the spectral index to be within 
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2(7 of the WMAP5 preferred value, = 0.963 ± 0.028 [16(.^ 

The motivation for measuring volumes which are consistent with inflation but not nec- 
essarily the particular experimental constraints we observe is the following. The volume of 
parameters satisfying some number of constraints would be a set of measure zero in the full 
parameter space. However we do not consider a model to be fine-tuned just because its 
parameters are fixed by certain measurements. Rather, it is the basic requirement of having 
a fiat enough potential to get at least 60 e-foldings of infiation (and not getting stuck in 
a local minimum preventing an exit from infiation) which underlies the apparent need for 
tuning that we are interested in. 

The actual objective function includes a few subtleties. For example, if a set of parameters 
does not give a monotonic potential nor satisfy the experimental constraints, then the relative 
volume can be defined to be zero. However, this gives no information to assist the Monte 
Carlo method in finding more favorable values since most points in the parameter space will 
have the same value of the objective function. In this case we therefore take the objective 
function (which is to be minimized) to be an empirical function of A^e, Vn, and n^: 

[ -On, On > 0. 

The empirical function is chosen in such a way as to help move the configuration from a 
nonmonotonic regime to a monotonic one, or from a regime not satisfying the experimental 
constraints to one which does. Once satisfactory parameters are found, then the relative 
volume is calculated, and maximizing 5n leads to parameter values which are less fine- 
tuned.^^ 

The Metropolis algorithm uses an artificial temperature T which randomly allows the 
objective function to sometimes increase rather than decrease; this is how it avoids getting 
stuck in a shallow local minimum rather than convergin g to some point closer to the global 



minimum. We used the downhill simplex method of Ref. [iTj as a generator of random steps, 
which moves the system's configuration by refiections, expansions, and contractions in an 
A^-dimensional simplex. As the temperature is lowered, the system relaxes to a minimum 



The power spectrum is insensitive to the location of the inflaton. For example, if one parametrization 
gives iVe = 55 and Us — 0.963 at V-r, = 2.41 x 10^^, then this parametrization (with a different inflaton 
location) will roughly give the same and Us at Vn = (2.41 ± 0.22) x 10"^ (WMAP5's 2a). Therefore 
we fix the normalization at WMAP's mean value for simplicity. 

The empirical function is not essential; it accelerates the search, but it has no effect on the results once a 
monotonic regime satisfying the constraints is found. The particular function used in this paper is given 
in Appendix iBl but one is free to design a different one. 
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which should be close to the global minimum. During this iterative process, a chain of 
accepted parameter values is generated, which allows one to make statistical statements 
about the probabilities of the parameters. 

B. Monte Carlo results 

Our starting point was the configuration given in the previous section, with the values 

Doi = 0.122 70, top = 15, Aq = 1, (p^ = 0.240 56, n = 8. (58) 

Since we had not yet varied enough parameters, this does not quite satisfy all the experimen- 
tal constraints: it has spectral index and normalization Ug = 0.949 and V-ji = 2.41 x 10~^ 
at = 78 e-foldings before the end of inflation, instead of at the COBE scale. Never- 
theless, since it is closer to the examples previously studied in the literature, we will use 
this as a reference point for comparison when assessing the improvement in fine-tuning. We 
also restarted the search using other initial conditions to avoid getting trapped in a local 
minimum. We accumulated approximately 200 chains containing more than 70 000 samples 
(each chain containing from 200 to 1000 samples). The control parameter T was decreased 
by a factor e after every m moves, over several orders of magnitude (for example, changing 
from 100 to 0.001) during the entire simulated annealing process. The computing time for 
each chain can be as fast as five hours (on a PC), or as slow as ten days, depending on the 
configuration of the annealing schedule and the initial conditions. 

To satisfy the experimental constraints we tried varying different combinations of the 
five independent parameters, {Doi,uf, Aq, (l)^,n). Generically one would expect that any 
two parameters would be uniquely fixed by the two constraints; thus to obtain chains, one 
should vary at least three at a time. Indeed, we found that the combinations {Dqi,ujp, Aq), 
{Dqi,u!p, 0^), or {Dqi, Up, n) were suitable for generating chains which evolved toward larger 
reduced volumes. In each case, the unvaried parameters take the values given in ( l58l) . which 
we refer to as the starting or fiducial point. The results are shown in Table [U where it can 
be seen that the reduced volumes of parameter space grow to values of order unity, starting 
from very small initial values. This indicates that there is essentially no fine-tuning at the 
optimal points. 

We can make more detailed statements about how the various parameters affect the 
power spectrum. For example, increasing up has the effect of reducing the overall scale of 
the potential, due to the e"^"^ dependence, while increasing Aq has the opposite effect. It is 
not surprising that the COBE normalization thus produces a strong correlation between Aq 
and LUp for the accepted parameter values, which can be seen in Fig. E] (left panel). Also 



We note that the first entry in Table U {S3 = 0.087) for {Doi,luf, Aq) is unphysical because we did not 
apply the constraint Aq < 1 (to avoid super-Planckian gaugino condensate scales) there. Starting from the 
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TABLE I: The relative volumes and reduced volumes for different combinations of the parameters, 
showing initial and final values along the chains. 



Parameters 


Starting point 


Optimal point 


{Doi,ujF,Ao) 

{Doi,UJF,(f)^) 

{Doi,iOF,n) 


9.0 X 10-'^ 0.96% 
5.2 X 10-1° 0.080% 
7.8 X 10-1° 0.092% 


0.087 44% 
0.027 30% 
0.020 27% 



shown in that Fig. [6] (right panel) is the correlation of 4D relative volumes with the value 
of Up for the Monte Carlo chains. The latter demonstrates that larger values of up are less 
likely than smaller ones, but only mildly so. 




FIG. 6: Left panel: scatter plot of accepted ujp and Aq values from search of 4D parameter space 
{Dqi,ujf, Aq, (p^). Right panel: correlation of ^4 and ljf- 

Furthermore, we find a degeneracy between 0^ and n on the shape of the potential: the 
effect of the discrete parameter n (the number of D7-branes in the stack) can be compen- 
sated by changing the value of 0^. Therefore the subset (-D01, ^^f, ^0, 0/^) gives an adequate 
representation of the possible potential shapes arising in the model. Fixing n = 8 to elim- 
inate the (p^-n degeneracy, the Metropolis algorithm finds a global minimum of /obj, hence 
a maximum of the relative volume ^4 = 0.052. Table HTl compares the optimal point to the 
fiducial point given in the previous section for the 4D search. We see that the fiducial point 

fiducial point (|58p . one cannot increase S3 by varying only {Dqi, ujf, Aq) while requiring j4o < 1. However 
the 4D parameter space search described below allowed us to find nontuned examples with ^0 < 1- 
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required tuning at the level of 0.5% per parameter; this number increases to 50% at the 
optimal point, so there is no more tuning. We also show the breakdown on a per-parameter 
basis: laij'pf' is the relative allowed width for the ith parameter. One might worry that the 
cruder statistic ^4 could hide severe tuning of some parameters by having very large values 
of 2ai/pf^ for others, to which inflation happened to be insensitive. However we see from 
Table HI] that this is not the case: the most sensitive parameter is Doi? which is still only 
tuned at the 20% level, in the region of the optimal parameter values. The last row of the 
table shows that the tuning per parameter is ameliorated by at least a factor of 100 (except 
for Aq, which did not require fine-tuning even at the starting point). 

TABLE II: Comparison of parameter values and degree of fine-tuning between the fiducial (starting) 
point and the optimal point of 4D parameter search. 



Configuration 


Dm 




Ao 








Fiducial point 


0.122 70 


15 


1 


0.240 56 


4.5 X 10" 


-10 0.46% 




1.3 X 10^3 


8 X 10"^ 


0.86 


5 X 10"^ 






Optimal point 


0.1976 


9.550 


0.007778 


0.5894 


0.052 


48% 




0.19 


0.52 


1.8 


0.3 






r / (0)1 


y 

opt 




140 

fid 


640 


2.1 


600 







Given that points in our chains tend to accumulate where the relative volume is bigger and 
the tuning problem is less severe, we can use the chains to define a probability distribution on 
the space of parameters, as well as on the predictions of the model for observable quantities. 
The correlation between relatively large volume and high density of models can be seen in 
the distribution of ^4, Fig. [71 The distributions for the spectral index (within the 2a range 
which we allowed around the WMAP5 central value) and the energy scale of inflation are 
shown in Fig. [HI We note there a preference for larger values of close to 0.99, while V^^^ 
is in the range (1.1-1.6) x lO'^Mpi = (2.7-3.9) x 10^^ GeV.^^ 



Recall that we impose three experimental constraints (Vtz, Ug, Ne) on the central point (pi ) of the volume, 
but we only require a monotonic potential for points in the volume, pf''^ — Ui < pi < pf^^ + cr^ . To solve 
the horizon problem, we need A'tot ~ 60; adding this requirement to the volume, the 2ai/p'f^ for 0^ in 
the table will be changed from 0.3 to 0.2. However, the modifications to the total volumes are small, 
5i = 0.036 and (5^ = 44%. And the conclusion remains true. 
Mpi is the reduced Planck mass, 2.44 x 10^^ GeV. 
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FIG. 7: Distribution of relative volume (^4 from the Monte Carlo chains. Notice the scale for 64 is 
logarithmic. 
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FIG. 8: Distribution of spectral index Ug and inflationary scale, based on Monte Carlo chains. 



C. An alternative measure 



Thus far, we have considered a well-defined and reasonable definition for the measure 
of fine-tuning on the parameter space: given a point which is consistent with the data, it 
measures the relative amount by which each parameter can be varied around this point and 
still be acceptable. However, one might argue that this hides fine-tuning problems in the 
case where the central values are orders of magnitude smaller than the largest theoretically 
allowed values (similar to the cosmological constant problem). In this section, we therefore 
consider an alternative measure, which is the absolute volume of parameter space consistent 
with measurements, Eq. fl5^ . instead of the relative volume. 

The absolute parameter volume as a measure of fine-tuning has some conceptual difficul- 
ties which are not present in the relative volume measure. First, it is possible that some 
parameters have a formally infinite range, such as the Kahler modulus up. One has to reg- 
ularize this divergence in order to make any sensible statements about probability. Second, 
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we must disentangle the issue of fixing parameter values through measurements from that 
of fine-tuning. For example, in the theory of quantum electrodynamics, we do not say that 
the charge of the electron is fine-tuned just because its value is experimentally known to a 
very high precision. The tuning problem we are concerned with is one where different free 
parameters have to be adjusted with respect to each other to a high precision in some artifi- 
cial way that had no a priori justification. Thus to implement the absolute volume measure, 
it makes sense to consider the volume only of some linear combinations of parameters which 
are not fixed by experimental constraints. This has the added advantage of possibly solving 
the first problem: we can use the experimental constraints to fix values of parameters which 
have no theoretical upper bound. 

There is a further shortcoming in the absolute volume measure when degeneracies exist 
between parameters. In the present case, the scale of the inflaton potential is determined by 
the product of two prefactors, and e"^"^^. As a result, a degeneracy between IuAq and 
ujp exists, as is confirmed by Fig. [6] (left panel). One can achieve a large absolute volume by 
trading up for Aq. For instance, we find models with super-Planckian values of Aq of order 
0(10^) {ujf ^ 30) having ctaq ~ O(IO^) and I/4 ~ 0{W),^^ whereas the absolute volume 
near the fiducial point is 2 x 10"^''; see Table UTTl But it would be misleading to claim that 
the model {ojf^Aq) ~ (30,10®) has less fine-tuning than that with {(jJf,Aq) ~ (10,0.01), 
because the large V4 is a result of the fact that the potential has different sensitivities to 
ujp and Aq. The relative volume, on the other hand, can balance the different sensitivities 
between different parameters. 

Despite these complications, it might still be of interest to know the ratio of the allowed 
volume to the total volume of the parameter space, 

A. = (59) 

i 

where Pi,max denotes the theoretical maximum of the parameter. As mentioned above, ujp 
has no maximum value, and the same is true of -Doi.max- Before dealing with this, let us 
consider just the numerator in fl59|) . Vat itself. We have repeated the Monte Carlo search 
of parameter space to find regions where Vat are maximized. The results are shown in 
Table IIIII There it can be seen that while Vn is only 2 x 10~^° for the fiducial (starting) 
point, it increases to 0.03 at the optimal point, which is different from the optimal point 
for the relative volume shown in Table HTl It is interesting that the improvement, 8 orders 
of magnitude, is the same as we obtained using the relative volume as the measure. Our 
conclusion that the degree of delicateness depends on the parameter values, and is greatly 
ameliorated in some regions of parameter space compared to others, thus holds for both the 
relative volume measure and the absolute volume measure. 
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^0 > 1 is theoretically disfavored; here we just use it to illustrate the effect of the degeneracy. 
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TABLE III: Comparison of parameter values and absolute allowed volume between the fiducial 
(starting) point and the optimal point of 4D parameter search. 



Configuration 


Dm 


UJF 


^0 






Fiducial point 


0.122 70 


15 


1 


0.240 56 


2 X 10"^° 


2ai 


1.6 X 10"^ 


0.012 


0.86 


1.2 X 10"^ 




Optimal point 


0.2186 


13.70 


0.6341 


0.4383 


0.03 




0.06 


11 


0.73 


0.063 




[o"4pt/Nfid 


380 


920 


0.85 


530 





Although the optimal case is much less fine-tuned than the fiducial point, we are still left 
with the question of how fine-tuned is it in an absolute sense, which would be answered if A^r 
was well defined. If we adopt the approach suggested above, of eliminating the parameters 
which have no maximum value {ujp and -Dqi), then we are left with A2 for the remaining 
parameters Aq and 0^. The maximum value of 0^ is 1, the maximum initial separation 
between the D3-branes and anti-D3-branes; similarly, Aq should be less than the Planck 
scale due to gaugino condensation. This gives A2 = 0.046 and ■\/A2 = 0.2 for the mean 
degree of tuning per parameter. This is comparable to the result we obtained using the 
relative volume measure. 

Alternatively, instead of eliminating the Dqi and ujp parameters, we could compare their 
widths to some "reasonable" or "generic" values. Since Dqi is the ratio of two brane tensions, 
Dqi ~ 1 would seem to be a generic value. As for up, even though in principle it could be 
arbitrarily large, in practice it is difficult to make it very large. It is a derived parameter, 
which depends upon Wq through Eq. f|T2l) . For example, if Wo,min ~ 10"''^, then cu^^max ~ 
100. We see that in fact it requires a high degree of fine-tuning of Wq to make ujp very large. 
Generously taking 100 to be its maximum value (even though this would be considered 
unreasonably large by most string theorists), we then find that A4 = 3 x 10~^, ~ 0.1, 
again a quite modest level of tuning per parameter. 
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VII. PROPERTIES OF THE OPTIMAL PARAMETER SET 

Having identified a favorable region in tfie space of the model parameters, we now 
consider a number of its detailed properties, including its consistency with string theoretic 
constraints, the shape of the inflationary trajectory, and sensitivity to initial conditions and 
the overshoot problem. We also explain a potential subtlety concerning the computation of 
the spectral index in the single-field approximation to the model. 

A. Microscopic parameter values 

The parameters {Doi,ujf, Ao,4'ii) which were convenient to vary in the potential are not 
all fundamental from the string theoretic point of view. We would like to determine the 
values of the microscopic stringy parameters which are compatible with the optimal point 
in Table [TTl which has 

^01 = 0.1976, oJir = 9.550, Ao = 0.007 778, 0^ = 0.5894, n = 8. (60) 

These yield the derived parameters uq = 9.644, s = 1.0094, and the corresponding observa- 
tional predictions are Vn = 2.41 x 10~^ for the primordial power spectrum and = 0.989 
for the spectral index, at = 52 e-foldings before the end of inflation, out of a total of 
A^tot = 134 e-foldings of inflation. It is straightforward to find a reasonable set of stringy 
parameters {N^, B4, Bq, Q^) giving the desired (0^, ujp)- For example, taking = 1.07 and 

= 10, we need B4 = 75.28, Bq = 1.0058. This set satisfies all the microscopic consistency 
conditions in Appendix lA 2[ 

One might however be concerned that such a small value of up as 9.6 is only marginally 
consistent with the need for control over higher derivative corrections to the low-energy 
effective action, which are suppressed by the large compactification volume (hence small 
curvatures). Figure [6] shows that larger values ofajp are indeed possible, up to a maximum 
of Up — 15; beyond this point, overly large values of Aq > 1 would have to compensate the 
reduction in the inflationary scale needed to get the right COBE normalization. Thus one 
can increase up to somewhat larger values, but at the expense of saturating the consistency 
condition Aq < 1, and somewhat increasing the degree of fine-tuning. However at up ~ 14.5, 
where the constraint Ao < 1 starts to become important, there is actually no fine-tuning: 
the reduced volume becomes ^4 ~ 0.021 at this point, as opposed to the optimal value 
^4 ~ 0.052. (Recall that 6'^ = = 38% quantifies the degree of tuning per parameter.) 



We use the optimal point given in Table HIl in this section, but the results also hold for the optimal point 
in Table HTH 
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B. Taylor expansion of DBI kinetic term 

Another point of consistency concerns the expansion of the DBI action for the inflaton 



kinetic term —a^T^yl — a~'^X'^ where T3 = m^/(87r^(yfs) is the 3-brane tension and a is the 
warp factor in the throat. Notice that the canonically normahzed inflaton field is ~ 
after Taylor expanding this expression. Estimating the inflation scale as l^^f^ ~ arris, we 
see that the criterion for being able to safely expand the DBI action into standard form 
is 0^ <^ T^a^ = V^nf/(8vr^(7s). On the other hand, the slow- roll equation of motion gives 
3H(j) ~ — V^. Since 3H^ ~ V^nf , this condition can thus be rewritten in terms of the slow-roll 
parameter e = as 

' « TTtV- (61) 

This is clearly satisfied in the present model, since as we will show, e ~ 10~^ at the horizon 
crossing. 



C. Initial conditions 

A potentially problematic aspect of the model is the need for special initial conditions, 
even if the potential itself is not finely tuned. Obviously, inflation takes place near the 
inflection point of the potential, so must not start lower than this point. But as was 
pointed out in Ref. [lo|, also should not start too much above the inflection point, because 
of the overshoot problem: the inflaton can gain so much speed that it quickly rolls past 
the inflection point without ever rolling slowly. For the optimal parameter set we consider, 
however, there is another consideration which prevents us from exploring the regime where 
overshoot would take place. This is because of the angular directions of the extra dimensions 
in the throat, which we have set to the values which minimize their potential. As shown in 
Ref. [sl, the positions of the angular minima flip when exceeds a certain critical value 0c, 
given by Eq. (C.24) of that paper. We would need to follow all the angular fields as well to 
investigate this regime quantitatively, which is beyond the scope of the present work. If the 
potential is assumed to have the same form for > 0c, we do observe overshooting, starting 
from initial conditions of order 100c, but since we do not trust the potential in this regime, 
no reliable statement about overshooting can be made in the present context. 

Nevertheless, we can quantify the range of initial conditions over which we get sufficient 
inflation and the potential is also valid: the allowed initial separation of the branes is 
X = 0.301 —>■ 0.665 for the fiducial point, while it is x = 0.029 —>■ 0.680 for the optimal 
point. We see that the allowed field range is expanded by a factor of 2 in the optimal 
case. Figure [9] shows the total number of e-foldings as a function of the initial conditions in 



these two cases. In an upcoming paper [l8|], we will give a more satisfactory solution to the 



problem of initial conditions in this model, based on the idea of Ref. 11|. 
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FIG. 9: The total number of e-foldings as a function of the initial condition Xj. The solid curve is 
for the fiducial parameters ([58]) . the dashed curve is for the optimal ones (j60|) . 



The present model has only one flat direction, so it is effectively a single-field inflation 
model. However, there is significant bending in the field space of uj and 0, as illustrated in 
Fig. [TOl In the left panel of Fig. [lOl the solid, smooth, leftmost curve is the "instantaneous 
minimum" uj^,{x), defined in Eq. (H9|) . while the wavy curve is the actual trajectory found by 
solving the equations of motion, given some initial displacement of the heavy field uo away 
from its instantaneous minimum. We used the initial condition uji = uj^,{xi) + 0.01, so there 
are oscillations at first which allow one to distinguish the two curves. The solid curve on 
the right is the effective single-field description given by Refs. 0,0]: 



As can be seen from the figure, it is not a very good approximation. Panda et al. [9] argued 
that using Eq. (162|) underestimates the total number of e-foldings by an order of magnitude. 
However, this does not invalidate the single-field description; the real instantaneous mini- 
mum uj^{x) does give a good approximation to the actual trajectory, once the oscillations 
have Hubble damped away. 

However a subtlety can arise in the computation of the spectral index when we use the 
slow-roll approximation ( l45l) . The slow-roll parameters along the adiabatic direction are 



related to those along the component field directions (f) = (f>^x and x = \/3/21nti; by Eqs. 



D. One field versus two: the spectral index 




(62) 
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FIG. 10: Left panel: actual inflationary trajectories in the valley-shaped potential, and a plot of 
approximate trajectory Eq. (j62p . Right panel: projection of the same trajectories onto the uj-x 
plane. 



(]A25P and (1A26P which we repeat here for convenience: 
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e^^e^ + e^, r]^^ ~ e^ + e^ i^^^*''' ^ ^''^^^ ^ "^^'^x) • (63) 

The problem occurs if we try to use these expressions to compute the spectral index Ug 
by assuming the inflaton rolls exactly along the instantaneous minimum u = uj^..{x). By 
definition, = e^^ = along this trajectory. However, if we neglect the terms proportional 
to in rj^^, we get a result which does not agree with computing the spectral index directly 
from d\nP/d\nk, Eq. (|40|). 

These conflicting results can be seen as the topmost (dot-dashed) and middle (dashed) 
curves, respectively, of Fig. [TT] (left panel), labeled as 77.^(0;*) and ns{Ckin)- The latter 
is based upon the approximation (|39|) for the power spectrum, and since it comes from 
directly differentiating Vn{k), it must be the correct result. The curve ns{uj^) significantly 
overestimates the spectral index as Ug = 1.065 near the inflection point {x = 0.029), whereas 
the correct value is = 0.989 for this example. Three approximations for the power 
spectrum itself, as a function of inflaton position x, are shown in the right panel of Fig. 
[TTl V-jiijP)) based on the slow-roll approximation (jSj), P7^(£kin) using (139!) along the exact 
trajectory, and Vniy^*) using (jHj) along the instantaneous minimum u = u!^:{x). It can be 
seen that they all agree quite well with each other near the inflection point, showing that 
u = Lj^{x) is indeed a good approximation. 

Nevertheless, the resolution of the problem is that the actual inflaton trajectory does 
not exactly follow the instantaneous minimum u = uj^:{x); it deviates slightly from this, 
like a race car on a banked curve. Therefore 7^ on the true trajectory (although it is 



As explained in Appendix |^ the slow-roll approximation is valid for both component fields since the 
inflaton itself is rolling slowly. 
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FIG. 11: Left panel: three approximations (described in the text) for the spectral index as a 
function of inflaton position x, indicating a problem with the slow-roll approximation on the 
instantaneous minimum trajectory, labeled ns(w*). Right panel: three approximations for the 
power spectrum as a function of x (see text). 



much smaller than e^), as plotted in the left panel of Fig. [121 and the two terms e^^xx 
V^V^rj^^/V"^ make an important contribution to rj^^. In fact, the right panel of Fig. [T2l shows 
that these two terms very nearly cancel e^r]^^ [l — Us = 0.011). The spectral index evaluated 
along the actual trajectory, but using the slow-roll formula f H5l) . is denoted ns{ip). Figure 
[TT] shows that it gives a good fit to the numerically computed index ?2<j(£kin) in the most 
important region, where most of inflation is taking place. Therefore the two-field slow-roll 
formula for Ug is a good approximation, but only if one uses the correct two-field trajectory 
and not the instantaneous minimum approximation. 
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FIG. 12: Left panel: the slow-roll parameters and e^-. is shown for both the exact trajectory 
and the approximation uj = a;*(2;). Notice that e^^ = in that approximation, so e^. is only shown 
for the exact trajectory. Right panel: the three terms which contribute to ij^^ in Eq. (j63p . 
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VIII. CONCLUSIONS 

In this paper we have made a detailed study of the warped D3-D3 inflation model, ac- 
counting for superpotential corrections, from a phenomenological perspective, but also with 
attention to the need for self-consistency from the string theoretical point of view. We ex- 
tended the model slightly by including uplifting from throats other than the inflationary one, 
which was necessary so that inflation could end with a nearly vanishing cosmological con- 
stant. We subsequently explored the parameter space using Monte Carlo methods, finding 
for the first time values which satisfy all theoretical and experimental constraints. Moreover 
we identified an optimal set of parameters in the vicinity of which there appears to be no 
need for fine-tuning. This arises in part because of a peculiar feature of the potential in this 
model: there are ranges of the parameter Dqi for which one finds a minimum number of 
e-foldings of inflation, because at the boundaries of such regions, a local minimum develops, 
giving rise to A^tot ^ oo. Toward this end, we defined the concept of a relative volume 
in an A^- dimensional parameter space consistent with successful inflation, and having the 
property that \/6n represents the average degree of fine-tuning of any parameter. We ob- 
tained an improvement by 8 orders of magnitude in ^4 for the 4-parameter subspace which 
was needed to fully describe the range of potential shapes in the model, relative to the 
fine-tuned example from the literature with which we started our Monte Carlo search. We 
conclude that this string theoretic inflation model is not as delicate as it first seemed in Ref. 

a. 

There are some caveats to this successful conclusion however. The value of the Kahler 
modulus at the optimal point, up = 9.6, is only marginally large enough to give one con- 
fidence that the low-energy effective field theory is not significantly perturbed by higher 
dimensional operators (coming from integrating out the extra dimensions) which are sup- 
posed to be suppressed by large up. On the other hand, since the scale of the potential goes 
like e~^'^^^o/-^Pb t^e COBE normalization demands that increasing up must be accompa- 
nied by an exponential increase in Aq. The latter is related to the energy scale A of gaugino 
condensation (or a Euclidean D3-brane) by A ~ A^, which should certainly not exceed the 
Planck scale, and likely should also lie below the warped string scale. This could be a source 
of theoretical tension for the model. 

On the phenomenological side, since we allowed for 2cr deviations of Ug from the WMAP5 
central value in our Monte Carlo search, we can detect a statistical preference of the model 
for larger values of Ug near 0.99. That is to say, the least fine-tuned models which we find 
correspond to such higher values of n^, and therefore if future data are shown to prefer lower 
values, it would be an indication disfavoring the model. In this way, our global search of the 
parameter space helps to provide an additional predictive tool which would not be available 
by simply finding a few sets of parameter values which were consistent with the data. This 
statement assumes that there is a landscape of string vacua which allows for nature to scan 
through the possible values of the parameters. 
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We also investigated an alternative measure of the degree of fine-tuning, i.e., the absolute 
allowed volume. While the motivation for this measure is different from that of the relative 
measure, the result is nearly the same — we found regions in the parameter space where 
the absolute allowed volume is increased by 8 orders of magnitude, and the overall degree 
of tuning per parameter is still at the 10% level. This finding reinforces our point that 
parameter regions exist where the severity of fine-tuning is greatly ameliorated. 

Note added. While we were finishing this work, Ref. 19|] appeared, where it was pointed 
out that generic deformations of the throat geometry can have a qualitatively similar effect 
to the superpotential corrections, in allowing for an inflection point in the potential. More 
recently Ref. 20|] appeared, which reaches a different conclusion than ours, finding too large 
a value of n,. Their work starts with a somewhat different model, based on the deformations 



discussed in Ref. 19j rather than superpotential corrections. But since the latter model is 
supposed to have qualitatively similar behavior to the former, the result of Ref. 20|] looks 
surprising. This discrepancy comes from the neglect of the Coulomb interaction in Ref. 20 1. 
We thank the authors of Ref. [sl for clarifying this point for us. 
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APPENDIX A: EQUATIONS 

Here we compile various formulas needed in the previous sections. 



1. LO<t:{x) and i^o 

The instantaneous minimum is defined through Eq. ( l49l) . This leads to the result 



1 



01 



Do 




- {2uf + 3) (^2 + cj, - ^wo^J 



3 / cx 



ng \9 



^uJo4>y + 6 = 0. 



Setting X = 0, we have the equation for ivq = ujq{s,ujf)'- 

2so;i.e2'^o-2<^^ - {2ujf + 3)(2 + uJo)e'^°-'^^ + 2cu^ + 7tUo + 6 = 0. 



(Al) 



(A2) 
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This equation generically has two solutions for ujq; the one closest to 00 p is the minimum 
while the one farther from cui? is at the maximum of the potential. 

If we apply the uplifting condition, V{0,u!o) = 0, then the parameter s is fixed and uq is 
a function of oup only; see Section IV C[ 

2. Microscopic constraints on the parameters 

There is a constraint on the allowed field range of the inflaton field jsl : 

If we consider < x < 1 only, then the parametrization Eq. (fT9l) requires > 1 and 
Qfj^ > 1. It is also required that S4 > 1. To satisfy the COBE normalization, one requires 
ujq < 0(30); otherwise the inflation scale will be too low. Using the fact that up — ujo leads 
to the further constraint N^/n < O(IO^). 



3. Potential and its derivatives 



We rewrite the potential as 

V = Vp + Vd = 



a|y4op 
3t/2 



where 



3 f cx 



LV~iJF^-l/n _i_ _^ _ ^3/2 



2cu + 6 - 2{2iup + 3)e'^~'^^^"^/" + — | x 



Vd = 2sujpe 



1 + 
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D- 



X^ 



aU 



Then we have 

dV ^ a\Ao\^ j ^Vf + Vd) 
duo 3t/2 

dV _ a|Aop 
'dx ~ 3[/2 

where 



2Vf + 2 



1 - {2iup + 3)e'^-'^^g 



3aU 



dx dx 



dx ng 
dVd 



V^{2up + 3)e--,-/" _3Vx/l cx\c 

9 \^ 9 J 9 



dx 



su;pe'''^''CDDiDi^x-° 1 + 



x^ 



(A4) 

(A5) 
(A6) 



, (A7) 
(A8) 



.,(A9) 
(AlO) 
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We need the second order derivatives when calculating the slow-roll parameters: 



where 



dx"^ 



8 / V dV\ 2a|AoP 
"aU [aU^du] ^ 3U^ 



{2Vf + e-2'^^2/n 



4u^(y^^^dV_2u^x 
Sail \ dx Sail 



4 + 3(2c<;^ + 3)e^-'^^^-^/"]}, 

(All) 
(A12) 
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(A13) 



(A14) 



and 



dxduj 



3aU 



' all ~^ du 



211 + -')? 

all J ox 



3f/2 \ dx 2ng 



2 - {2iUF + 3)e'^''^^g 



uj—ujp —1/n 



(A15) 



4. Power spectrum 



The kinetic part of the Lagrangian is 

a2[/2 /I 



kin 



-TT, 



2uJoUJ(f)l 



(A16) 



Its derivative is given by 



(i£kin 4£kin / dto ^o(PuX dx\ ^ a^U"^ 



dN 



aU \dN 



3 dN 



4 



271 dn^ 



+ 



dTTx TTx dui 



3 dN ujQUj(t)l \ dN 2uj dN 



(A17) 



where dx/dN, du/dN, dTT^/dN , and dTc^^/dN are obtained by solving the equations of mo- 
tion. From Eq. fl28|) . we have 



d\nH 

dN 



dCv\ 



kin 



dx ^duj\ 

, K h V,, — . 

\ dN ""dN ""dN 



(A18) 



The spectral index can be calculated through Eq. (j40 
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5. Slow-roll parameters 



The slow-roll parameters for a generic field (f) are defined as 

1 fVs^ ' 



= 2 ) ' (A19) 



V 



(A20) 



where V^^ = d^V/dcj?. Since [15 1 



= 1/0cos^ + y^sin^, (A21) 
V^i/i = V^iji COS ^ + V^^ sin^ ^ + 1/^^ sin 29, (A22) 

we have 

e, = + + (A23) 

+ X 

In general, one should not apply the slow-roll approximation to all the terms in these ex- 
pressions since one linear combination of the fields is heavy. However, we are interested in 
trajectories along the flat direction of the potential, after any oscillations in the steep di- 
rections have Hubble damped away. In this case, the slow-roll approximation is valid along 
both of the field components, and we can write 

ev, ^ e^ + ^x, (A25) 
^"^"^ - e^ + e^ V ^^'^'^ ^ ^^^^^ ^ 'yr^'t'x ) ■ (A26) 
Using the relation between (x, a;) and (0, x); we have 

= 4ex, (A27) 
0M 

ex = (A28) 

= 4'7xx, (A29) 



^xx = 3^ (y + ^'^^^ ' ^^'^^^ 

V^x = \j^j-Jlxu- (A31) 

The fact that the slow-roll approximation for the spectral index agrees with the numerical 
computation of Ug provides further evidence for the validity of the approximation. 



34 



APPENDIX B: OBJECTIVE FUNCTION FOR MONTE CARLO METHOD 

In Section |VT] we described the Metropolis algorithm for finding parameters which maxi- 
mize the relative volume. The strategy is to minimize an objective function. A naive choice 
of the objection could be the negative relative volume: 

/obj = -Sn. (B1) 

However, the concept of the relative volume is only valid when the potential is monotonic, 
so the objective function becomes zero in the nonmonotonic regions, 

jo, V nonmonotonic, , , 

-On, On > 0. 

But this is not a useful choice, since most points in the iV-dimensional space give non- 
monotonic potentials. Supposing that we start with an A^-dimensional simplex, if we do 
not choose the + 1 initial vertices carefully, then there is a good chance that all + 1 
vertices correspond to nonmonotonic potentials. Therefore all the initial vertices have the 
same value of the objective function, and it may take a long time for the code to escape 
from such a region, since, apparently, there is no downhill direction. This problem might 
be solved by selecting good initial vertices; however, it is not practical because, before first 
getting a few successful chains, one does not know the configuration of the A^- dimensional 
parameter space. 

To avoid this situation, we need another objective function for the nonmonotonic regime. 
The primary characteristic of a nonmonotonic potential, having a local minimum, is that 
the number of e-foldings diverges as the infiaton gets stuck in the minimum. However, in 
the Runge-Kutta method (with adaptive stepsize control) which we use to solve the infiaton 
equations of motion, a parameter hmin controls the minimum stepsize, which causes the 
evolution to end after a finite number of e-foldings. This number will tend to be larger 
for a shallow local minimum than for a deep one. Thus, to assist the program in finding 
parameters that move away from a local minimum, a good choice for the objective function 
in the nonmonotonic regime is 

\ N' - A^tot, V nonmonotonic, 

[ -On, On > U. 

where A^' is some large number which is not supposed to be attained in the code. 

The above objective function should work. However, we not only want the central point to 
give a monotonic potential, but we also need it to satisfy the experimental cosmic microwave 
background constraints. The next step is to find the value of the infiaton field x where the 
COBE normalization is satisfied 'P7^(xcoBE) = 2.41 x 10^^. (Except for the case whose 
whole power spectrum is less than the COBE normalization scale, this point xcobe can 
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always be found.) To make this correspond to the right scale of wave numbers, we also need 
the number of e-foldings at this field value to satisfy 50 < A^cobe < 60. This leads to the 
further refinement: 



/obj 



50 - A^COBE, 

-/VcoBE - 60, 



V nonmonotonic, 
A^COBE < 50, 
A^COBE > 60, 

Sn > 0. 



(B4) 



If the constraint on Nqq^e is satisfied, then the next step is the constraint on the spectral 
index. We require it to be within 2a of WMAP5's mean value, so we add 



/c 



A^' 
50- 



obj 



■ A^tot, 

A'COBE, 
A^COBE - 60, 
|n, -0.963|, 

-Sn: 



V nonmonotonic, 
iVcoBE < 50, 

A^COBE > 60, 

50 < iVcoBE < 60, 

0.935 < n, < 0.991. 



(B5) 



Finally, we must avoid overlapping conditions between the above cases. We thus modify 
it to 

' max(50,Ar'- ATtot), 
max(2,50 - A^cobe), 
min(50,max(2, A^COBE - 60)), 
\ns - 0.963|, 



/obj 



JN, 



V nonmonotonic, 
A'^COBE < 50, 
iVcoBE > 60, 
50 < A^coBE < 60, 
0.935 < < 0.991, 



(B6) 



where the number 50 is the minimum number of e-folding we need, while the number 2 is 
somewhat arbitrary (we suppose that \ns — 0.963| < 2). 
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